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ABSTRACT 

We investigate various galaxy occupation statistics of dark matter halos using a large galaxy group 
catalogue constructed from the Sloan Digital Sky Survey Data Release 4 (SDSS DR4) with an adaptive 
halo-based group finder. The conditional luminosity function (CLF), which describes the luminosity 
distribution of galaxies in halos of a given mass, is measured separately for all, red and blue galaxies, as 
well as in terms of central and satellite galaxies. The CLFs for central and satellite galaxies can be well 
modelled with a log-normal distribution and a modified Schechter form, respectively. About 85% of 
the central galaxies and about 80% of the satellite galaxies in halos with masses Mh > lO^^/i~^M0 are 
red galaxies. These numbers decrease to 50% and 40%, respectively, in halos with Mh ~ 10^^ /i~^Mq. 
For halos of a given mass, the distribution of the luminosities of central galaxies, L^, has a dispersion 
of about 0.15 dex. The mean luminosity (stellar mass) of the central galaxies scales with halo mass 
as Lc oc Mjl'^"^ (M,^c oc M^-^^) for halos with masses M > 10^^-^ H'^Mq, and both relations are 
significantly steeper for less massive halos. We also measure the luminosity (stellar mass) gap between 
the first and second brightest (most massive) member galaxies, logLi — logL2 (logM<,_i — logM*_2)- 
These gap statistics, especially in halos with Mh ^ 10^^ '^ h~^MQ, indicate that the luminosities of 
central galaxies are clearly distinct from those of their satellites. The fraction of fossil groups, defined 
as those groups with logLi — \0gL2 > 0.8, ranges from ^ 2.5% for groups with Mh ~ 10^''' H^^Mq to 
18-60% for groups with Mh ~ 10^'^ H^^Mq. The number distribution of satellite galaxies in groups 
of a given mass follows a Poisson distribution, in agreement with the occupation statistics of dark 
matter sub-halos. This provides strong support for the standard lore that satellite galaxies reside in 
sub-halos. Finally, we measure the fraction of satellites, which changes from ~ 5.0% for galaxies with 
°-^Mr - 5 log /i 22.0 to - 40% for galaxies with " ^Mr - 5 log /i 17.0. 

Subject headings: dark matter - large-scale structure of the universe - galaxies: halos - methods: 
statistical 



1. INTRODUCTION 

In recent years, the halo occupation distribution and 
conditional luminosity function have become powerful 
statistical measures to probe the link between galaxies 
and their hosting dark matter halos. Although these 
statistical measures themselves do not give physical ex- 
planations of how galaxies form and evolve, they pro- 
vide important constraints on various physical processes 
that govern the formation and evolution of galaxies, such 
as gravitational instability, gas cooling, star formation, 
merging, tidal stripping and heating, and a variety of 
feedback processes. In particular, they constrain how 
their efficiencies scale with halo mass. 

The halo occupation distribution (hereafter HOD), 
P{N\M), which gives the probability of finding N galax- 
ies (with some specified properties) in a halo of mass M, 
has been extensively used to study the galaxy distribu- 
tion in dark matter halos and galaxy clustering on large 
scales (e.g. Jing, Mo & Borner 1998; Peacock & Smith 
2000; Seljak 2000; Scoccimarro et al. 2001; Jing, Borner 
& Suto 2002; Berlind & Weinberg 2002; Bullock, Wech- 
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sler & Somerville 2002; Scranton 2002; Kang et al. 2002; 
Marinoni & Hudson 2002; Zheng et al. 2002; Maghoc- 
chetti & Porciani 2003; Berlind et al. 2003; Zehavi et 
al. 2004, 2005; Zheng et al. 2005; Tinker et al. 2005). 
The conditional luminosity function (CLF), ^{L\M)dL, 
which refines the HOD statistic by considering the aver- 
age number of galaxies with luminosity L ± dL/2 that 
reside in a halo of mass M, has also been extensively 
investigated (Yang, Mo & van den Bosch 2003; van den 
Bosch, Yang & Mo 2003; Vale & Ostriker 2004, 2006; 
Cooray 2006; van den Bosch et al. 2007a) and has been 
applied to various redshift surveys, such as the 2dFGRS, 
the Sloan Digital Sky Survey (SDSS) and DEEP2 (e.g. 
Yan, Madgwick & White 2003; Yang et al. 2004; Mo et al. 
2004; Wang et al. 2004; Zehavi et al. 2005; Yan, White 
& Coil 2004). These investigations demonstrate that the 
halo occupation statistics are very useful in establishing 
and describing the connection between galaxies and dark 
matter halos. Furthermore, they also indicate that the 
galaxy/dark halo connection can provide important con- 
straints on cosmology (e.g., van den Bosch, Mo & Yang 
2003; Zheng & Weinberg 2007). Finally, the HOD/CLF 
framework also allows one to split the galaxy population 
in centrals and satellites, and to describe their properties 
separately (e.g. Cooray 2005; White et al. 2007; Zheng 
et al. 2007). 

As has been pointed out in Yang et al. (2005c; here- 
after Y05c), a shortcoming of the HOD/CLF models is 
that the results are not completely model independent. 
Typically, assumptions have to be made regarding the 
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functional form of either P{N\M) or ^{L\M). More- 
over, in all HOD/CLF studies to date, the occupation 
distributions have been determined in an indirect way: 
the free parameters of the assumed functional form are 
constrained using statistical data on the abundance and 
clustering properties of the galaxy population. One may 
hope to circumvent this shortcoming by directly mea- 
sure the dark matter distribution around galaxies. Such 
measurements can in principle be obtained through grav- 
itational lensing and X-ray observations. However, both 
methods are hampered by requirements on the data qual- 
ity and uncertainties related to the interpretation of the 
data. For instance, weak lensing measurements, which 
requires high-quality imaging, typically needs to resort 
to the stacking of many lens galaxies in order to get a 
detectable signal, but this stacking severely complicates 
the interpretation in terms of the halo masses of the lens 
galaxies. In the case of X-ray observations, robust con- 
straints can only be obtained for massive clusters, but 
even here the interpretation of the data can be compli- 
cated due to the presence of substructure and deviations 
from hydrostatic equilibrium. An alternative method to 
directly probe the galaxy - dark halo connection is to use 
galaxy groups as a representation of dark matter halos 
and to study how the galaxy population changes with 
the properties of the groups (e.g., Y05c; Zandivarez et 
al. 2006; Robotham et al. 2006; Hansen et al. 2007). 

Recently, we have constructed a large galaxy group 
catalogue based on the Sloan Digital Sky Survey Data 
Release 4 (SDSS DR4), using an adaptive halo-based 
group finder (Yang et al. 2007; Paper I; Y07 hereafter). 
Detailed tests with mock galaxy catalogues have shown 
that this group finder is very successful in associating 
galaxies according to their common dark matter halos. 
In particular, the group finder performs reliably not only 
for rich systems, but also for poor systems, including iso- 
lated central galaxies in low mass halos. This makes it 
possible to study the galaxy-halo connection for systems 
covering a large dynamic range in masses. Various obser- 
vational selection effects have been taken into account, 
especially the survey edge effects and fiber collisions. The 
halo masses for the groups are estimated according to 
the abundance match, using the characteristic group lu- 
minosity and stellar masses (see fj2] below). According 
to tests with mock galaxy catalogues, the halo masses 
are estimated with a standard deviation of about 0.3 
dex. With these well-defined galaxy group catalogues, 
one can not only study the properties of galaxies in dif- 
ferent groups (e.g. Y05c; Yang et al. 2005d; CoUister 
& Lahav 2005; van den Bosch et al. 2005; Robotham et 
al. 2006; Zandivarez et al. 2006; Weinmann et al. 2006a, b; 
van den Bosch et al. 2007b; Mcintosh et al. 2007), but 
also probe how dark matter halos trace the large-scale 
structure of the Universe (e.g. Yang et al. 2005b, 2006; 
Coil et al. 2006; Berlind et al. 2007; Wang et al. 2007 
in preparation). In this paper, which is the second in a 
series, we use the SDSS DR4 group catalogue to probe 
various occupation statistics and measure the CLFs for 
different populations of galaxies. 

This paper is organized as follows: In Section [2] we 
describe the data (galaxy and group catalogues) used 
in this paper. Section [3] presents our measurement of 
the CLFs for all, red and blue galaxies. Sections [H [5] 
and [6] describe the properties of central galaxies, the 
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Fig. 1. — The color-magnitude relation for galaxies in our group 
sample. The open circles indicate the Gaussian peaks of the bi- 
normal distribution of galaxies in each luminosity bin. The solid 
dots indicate the corresponding averages of the two Gaussian peaks. 
The solid line is the best-fit quadratic relation to these averages 
(see eq. [l]), which we use to split the galaxies into red and blue 
population (color-coded accordingly). 

halo occupation statistics and the fraction of satellite 
galaxies, respectively. Finally, we summarize our results 
in Section [T] Throughout this paper, we use a ACDM 
'concordance' cosmology whose parameters are consis- 
tent with the three-year data release of the WMAP mis- 
sion: rim = 0.238, f^A = 0.762, n, = 0.951, h = 0.73 and 
as = 0.75 (Spergel et al. 2007). 

2. DATA 

The analysis presented in this paper is based on the 
SDSS DR4 galaxy group catalogue constructed by Y07 
using an adaptive halo-based group finder (Yang et 
al. 2005a). The related galaxy catalogue is the New 
York University Value-Added Galaxy Catalogue (NYU- 
VAGC; Blanton et al. 2005b), which is based on the 
SDSS Data Release 4 (Adelman-McCarthy et al. 2006), 
but with an independent set of significantly improved 
reductions. From this catalogue we select all galax- 
ies in the Main Galaxy Sample with redshifts in the 
range 0.01 < z < 0.20 and with a redshift complete- 
ness C > 0.7. As described in Y07, three group samples 
were constructed from the corresponding galaxy samples: 
Sample I, which only uses the 362356 galaxies with mea- 
sured r-band magnitudes and redshifts from the SDSS, 
Sample II which also includes 7091 galaxies with SDSS 
r-band magnitudes but redshifts taken from alternative 
surveys, and Sample HI which includes an additional 
38672 galaxies that lack redshifts due to fiber collisions 
but that are assigned the redshifts of their nearest neigh- 
bors (cf. Zehavi et al. 2002). Unless stated otherwise, 
our analysis is based on Sample II. For comparison, we 
also present some results obtained from Sample III. 

The magnitudes and colors of all galaxies are based on 
the standard SDSS Petrosian technique (Petrosian 1976; 
Strauss et al. 2002), have been corrected for galactic ex- 
tinction (Schlegel, Finkbeiner & Davis 1998), and have 
been if-corrected and evolution corrected to z = 0.1, 
using the method described in Blanton et al. (2003a, b). 
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Fig. 2. — The conditional luminosity functions (CLFs) of galaxies in groups of different mass bins. Symbols correspond to the CLFs 
obtained using A/^ as halo mass (estimated according to the ranking of the characteristic group luminosity), with solid and open circles 
indicating the contributions from central and satellite galaxies, respectively. The errorbars reflect the 1-cr scatter obtained from 200 
bootstrap samples. The solid lines indicate the related best-fit parameterizations using equation [2]. For comparison, we also show, with 
dashed lines, the CLFs obtained using Ms as halo mass (estimated according to the ranking of the group's characteristic stellar mass). 
Results shown in this plot are obtained from Sample II. 



We use the notation ^-^Alr to indicate the resulting abso- 
lute magnitude in the r-band. Stellar masses, indicated 
by M^, for all galaxies are computed using the relations 
between stellar mass-to-light ratio and — r) color of 
Bell et al. (2003; see Y07 for details). 

In this study we separate galaxies into red and blue 
subsamples according to their bi-normal distribution in 
the °-^{g - r) color (Baldry et al. 2004; Blanton et 
al. 2G05a; Li et al. 2006). Fig. [1] shows the color- 
magnitude distribution of the galaxies in our Sample II 
(dots) together with the two peak values of the bi-normal 
distribution in each absolute magnitude bin (open cir- 
cles) (Cheng Li; private communication). The galaxies 
are separated into red and blue subsamples using the 
solid line, which is the best fit to the average of the two 
peak values in each absolute magnitude bin: 

-r) = L022 - 0.0651x - O.OOSllx^ , (1) 



where x= '^ iM^ - 5 log/i 23.0. 

For each group in our catalogue we have two estimates 
of its dark matter halo mass M^. (1) M^, which is based 
on the ranking of the characteristic group luminosity 
Li9.5 , and (2) M5, which is based on the ranking of the 
characteristic group stellar mass Mstciiar, respectively^. 
As shown in Y07, these two halo masses agree reason- 
ably well with each other, with a scatter that decreases 
from ^0.1 dex at the low-mass end to ~ 0.05 dex at the 
massive end. Detailed tests using mock galaxy redshift 
surveys have demonstrated that the group masses thus 
estimated can recover the true halo masses with a l-cr 
deviation of 0.3 dex, and are more reliable than those 
based on the velocity dispersion of group members (Y05c; 

^ -L19.5 and MgtgUar S'l's, respectively, the total luminosity and 
total stellar mass of all group members with "^'^Mr — 51og/i < 
-19.5. 
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Fig. 3. — Similar to Fig. [2] but here we show the CLFs for red (dashed Unes) and blue (dotted lines) galaxies, for groups with halo 
masses Mi^. The solid lines indicate the best-fit parameterizations for the CLFs of red galaxies. In both cases the central and satellite 
components of the CLFs are indicated separately. For clarity, the errorbars, again obtained using 200 bootstrap samples, are only shown 
for the red galaxies. 



Weinmann et al. 2006; Berlind et al. 2006; Y07). Note 
also that survey edge effects have been taken into account 
in our group catalogue: groups that suffer severely from 
edge effects (about 1.6% of the total) have been removed 
from the catalogue. In most cases, we take the brightest 
galaxy in the group as the central galaxy (BCG) and all 
others as satellite galaxies. In addition, we also consid- 
ered a case in which the most massive galaxy (in terms 
of stellar mass) in a group is considered as the central 
galaxy (MCG). Tests have shown that for most of what 
follows, these two definitions yield indistinguishable re- 
sults. Whenever the two definitions lead to significant 
differences, we present results for both. 

3. THE CONDITIONAL LUMINOSITY FUNCTIONS 
FOR ALL, RED AND BLUE GALAXIES 

The conditional luminosity function (CLE) of galax- 
ies in dark halos, $(i|A/), which describes the average 
number of galaxies as a function of galaxy luminosity in 



a dark matter halo of a given mass, plays an important 
role in our understanding of how galaxies form in dark 
matter halos (Yang, Mo & van den Bosch 2003; van den 
Bosch, Yang & Mo; 2003; Cooray 2006; Vale & Ostrikcr 
2006; van den Bosch et al. 2007a). We now use our group 
catalogue to directly determine $(L|M). 

The CLE can be estimated by directly counting galax- 
ies in groups. Eor a given galaxy luminosity L, there 
is a limiting redshift, zl, beyond which galaxies with 
such a luminosity are not included in the sample. As 
discussed in Y07, the group catalogue is complete to a 
certain limiting redshift for a given halo mass. In order 
to estimate the CLE, ^{L\M) at a given L, we only use 
groups that are complete to the redshift limit z^. The 
CLE is obtained by simply counting the average num- 
ber of galaxies (in luminosity bins) in groups of a given 
M. We show in Eig. [2] the resulting CLEs for groups 
of different masses. The contributions of central and 
satellite galaxies are plotted separately. Eor compari- 
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Fig. 4. — The best fit parameters {(p* upper row, a* second row, Lc third row, and <Tc bottom row) to the CLFs shown in Figs. [2] and 
[S] as functions of halo mass. Panels on the left, in the middle, and on the right show results for all, red, and blue galaxies, respectively. 
Since we have two different halo mass estimators (M^ and A/g) and two main group samples (II and III), we have obtained CLFs for four 
different combinations of sample and group mass estimator. The results for all four combinations are shown using different symbols and 
line-styles, as indicated. The errorbars in the first two and last rows indicate the l-tr variances obtained from our 200 bootstrap samples. 
In the third row of panels, however, the errorbars correspond to the log-normal scatter, CTc, shown in the bottom row of panels. For clarity 
the errorbars are only shown for the 'Af^-Sample IF case, but they are very similar for the other four cases. 



son, results obtained using both and Ms are shown 
as symbols and dashed lines, respectively. Overall, these 
two halo masses give consistent results, except that the 
Mi-based CLF of the central galaxies in low mass ha- 
los is more peaked than the Mg-based CLF (see the 
lower right-hand panel). The errorbars shown in each 
panel correspond to 1-a scatter obtained from 200 boot- 
strap samples of our group catalogue. The CLF for the 
total population (centrals plus satellites) matches rea- 
sonably well the Schechter form down to halo masses of 
M ~ 10^^'^ H'^Mq. For less massive halos, however, 
there is a prominent peak in the CLF at the bright end 
due to the contribution of central galaxies, which makes 
the total CLF deviate significantly from the Schechter 
form. As discussed in Y05c, in low-mass halos where the 
group characteristic luminosity (iig.s) is dominated by 
the central galaxy, the Liq^^-Ml conversion can produce 
an artificial peak in the CLF at the bright end. However, 
comparing the results obtained here with the test results 
shown in Fig. 9 of Y05c indicates that the strong peak 
in the lower right panel of Fig. [5] cannot be entirely due 
to the Li9.5-Ml conversion. Therefore, in what follows, 



we will model the CLF for central and satellite galaxies 
separately. 

Fig.[3]shows the CLFs separately for red (dashed lines) 
and blue (dotted lines) galaxies. Note that massive halos 
clearly contain more red galaxies than blue galaxies (both 
centrals and satellites), while the opposite applies to low 
mass halos. The overall CLF shapes for red and blue 
galaxies, however, are remarkably similar. 

To quantify the CLFs, we fit each of them with the 
following model. We write the total CLF as the sum of 
the CLFs of central and satellite galaxies: 



$(L|M) = $ccn(i|Af) + $sat(i|M) . 



(2) 



We assume the contribution from the central galaxies to 
be a lognormal: 



$ccn(i|M) = 



A 



2n(T, 



■exp 



(logi-logL,)2 



2^2 



(3) 



where A is the number of central galaxies per halo. Thus, 
A = 1 for all galaxies, A = /rod (with /red being the 
fraction of red central galaxies) for red galaxies, and 
A ~ 1 — /red for blue galaxies. Note that log Lc is, by 
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Fig. 5. — Fraction of red galaxies among central (left-hand panel) and satellite (right-hand panel) galaxies. Solid and dashed lines show 
results for groups with halo mass estimators and Ms, respectively. For comparison, we show also the fractions of red central galaxies 
obtained from the CLFs (Eq. [Sj as open circles (for Mi^) and solid squares (for Mg). Errorbars (shown only for -A/^) are obtained from 
the 200 bootstrap samples. See text for a detailed discussion. 



definition, the expectation value for the (10-based) loga- 
rithm of the luminosity of the central galaxy; 



logic 



$ccn(i|Af)logidL, 



(4) 



and that CTc = a{\ogLc). For the contribution from the 
satellite galaxies we adopt a modified Schechter function: 



$.at(i|A/)=4* ( 



exp 



(5) 



Note that this function decreases faster at the bright end 
than a Schechter function. This CLF parameterization 



has a total of five free parameters: i^, 



a* and L* 



In practice, we find that log Lf. ~ log L*-|-0.25 to good ap- 
proximation. We therefore adopt logLc = fog^s + 0.25, 
throughout, which reduces the number of free parameters 
to four. 

For all the CLFs measured above, the best fit parame- 
ters, using the model described by equation[2l are shown 
in Fig [31 Results are shown separately for all (left col- 
umn panels), red (middle column panels) and blue (right 
column panels) galaxies. The error-bars in the upper two 
rows reflect the l-cr scatter obtained from our 200 boot- 
strap samples. Here we also compare results obtained 
using Samples II and III. Note that Sample II does not 
include any galaxies missed due to fiber collisions, while 
Sample III includes all such galaxies by assigning each of 
them the redshift of its nearest neighbor. Although this 
fiber collision correction works well in roughly 60 percent 
of all cases, the remaining 40 percent are assigned red- 
shifts that can be very different from their true values 
(Zehavi et al. 2002). Samples II and III should therefore 
be considered as two extremes as far as a treatment of 
fiber-collisions is concerned. Given that we also have two 
kinds of halo masses. Ml and Mg, we have a total of four 
different combinations for which we have determined the 
CLF. The results for all four cases are shown as different 
symbols in each of the panels of Fig ID As an illustra- 
tion of how well the model fits the data, the solid lines in 
Figs. [2] and [3] indicate the corresponding best-fit models. 



The upper row of Fig. |4] shows the best fit normal- 
ization of the CLF for satellite galaxies, which describes 
the average number of satellite galaxies with luminosity 
^ L* in a group of a given halo mass. As expected. Sam- 
ple III gives a higher 0*, especially for low-mass groups. 
Comparing (f>* for red (upper middle panel) and blue (up- 
per right panel) galaxies, one sees that the fraction of 
red satellites increases with halo mass. The second row 
shows the faint end slopes of the CLFs, a*. In massive 
haloswithMh ^ 10^^/i~^Mq, a* decreases (i.e., becomes 
more negative) with increasing halo mass, both for red 
and for blue galaxies. In halos with Mh < 10^^ ft.~^M0, 
however, a* decreases with increasing halo mass for red 
satellites, but increases with Mh for blue satellites, while 
the faint-end slope for the entire satellite population (red 
and blue combined) is roughly constant at ~ —1.1. The 
third row of panels shows that log Lc increases with halo 
mass, for both red and blue centrals. A more detailed 
discussion regarding this and other properties of central 
galaxies will be presented in Section |4]) . Finally, the last 
row of Fig 0] shows the width of the log-normal CLF of 
central galaxies. For the combined sample of red and blue 
galaxies we find an average value of CTc = o'(log L^) ~ 0.15 
This is in good agreement with constraints obtained by 
Yang et al. (2003) and Cooray (2006) from clustering 
and abundances of galaxies in the 2dFGRS and SDSS, 
respectively. However, Zheng et al. (2007) found, based 
on HODs models for the SDSS and DEEP2, that the 
log-normal width increases from 0.13 for massive ha- 
los with Mh ^ 10^^-^ h~^MQ to ~ 0.3 for low mass halos 
with Mh ~ 10"-^ /i-^Mq. This is quite different from 
our results, which suggest that a(\ogLc) decreases with 
decreasing halo mass. However, this is most likely an 
artefact of the method that we used to assign halo masses 
to our groups. As mentioned above, our halo masses are 
based on the ranking of either L19.5 or Mstdiar, which 
implies that we have assumed a one-to-one relation be- 
tween halo mass and these two indicators. At the low 
mass end, both L19.5 and Mstciiar are expected to be 
strongly correlated with the luminosity of corresponding 
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central galaxies. Therefore the values of CTc obtained here 
should be considered as lower limits, in particular for low 
mass halos. 

For reference, Table [T] lists the average values of the 
CLF fitting parameters obtained from the four combi- 
nations of Samples II and III and group masses and 
Ms- The errorbars indicate the scatter among these four 
combinations or the scatter obtained from the 200 boot- 
strap samples, whichever is larger. 

As shown in Weinmann et al. (2006b) and Baldry et 
al. (2006), the red and blue fractions of galaxies as func- 
tion of halo mass provide important constraints for mod- 
els of galaxy formation. Fig. [5] shows the red fractions 
of centrals (left-hand panel) and satellites (right-hand 
panel) as function of halo mass obtained from our group 
catalogue (solid and dashed lines). These fractions are 
obtained by simply dividing the number of red centrals 
(satellites) by the number of all centrals (satellites) in 
a given halo mass bin. For both centrals and satellites 
the red fraction increases with increasing halo mass. For 
comparison, results are shown for both (solid lines) 
and Ms (dashed lines). These two different halo mass 
estimates give quite different results for central galaxies, 
especially for halos with Mh ~ 10^^ h~^MQ. The origin 
of this discrepancy can be understood in terms of Fig. [TJ 
In this color-magnitude diagram, the galaxies are sepa- 
rated into red and blue populations. If we would convert 
this color-magnitude diagram into a color-stellar mass di- 
agram, the red and blue populations would shift slightly 
towards higher and lower stellar masses, respectively (ac- 
cording to the Eq. 2 in Y07). Therefore, a certain upper 
percentile of galaxies that is ranked according to stellar 
mass will have a larger red fraction than the same per- 
centile ranked according to luminosity. Because of the 
tight correlation between the luminosity of the central 
galaxy and Ml , and between the stellar mass of the cen- 
tral galaxy and Ms, it is clear that centrals in a given 
bin of Ms have a larger red fraction than those in the 
same bin of M^. 

We can also determine the red fractions of centrals 
from our best-fit CLF parameterizations (equation [3]). 
As shown above, the parameterization of ^ccn{L\-^^) 
volves the parameter /red, which describes the fraction 
of central galaxies that are red. For comparison, the 
symbols in the left-hand panel of Fig. [5] show the best- 
fit values of /red obtained from the CLFs for the 
masses (open circles) and Ms masses (solid squares). It 
is reassuring that these best fit results match the direct 
determination of the red fractions reasonably well. The 
small differences between the direct measurement (lines) 
and CLF measurement (symbols) owe to the fact that 
in the latter case we have normalized the galaxies using 
the number of groups within the redshift limit (see 
Section [3]). 

The results shown in the left panel indicate that more 
than 80% of the central galaxies in halos more massive 
than ~ 10^'' h~^MQ are red galaxies, while in smaller 
halos with masses ^ 10^^ h~^MQ, less than 50% of the 
centrals are red. As discussed in Y07, AIs may represent 
the true halo mass better than M^. The result based 
on Ms shows that /red decreases rapidly with decreas- 
ing halo mass at Afh ^ 10^"^ h~^MQ. The right panel 
of Fig. [5] shows the fraction of red satellite galaxies as 
a function of halo mass. As one can see, for satellite 



galaxies /red increase steadily from about 40% for halos 
with Mh ~ 10^^ h~^MQ to about 80% in massive halos 
with Mh ~ 10^5 /i-^Mq. In van den Bosch et al. (2007b) 
we have used this information to constrain the efficiency 
with which the star formation of galaxies is quenched 
once they become a satellite galaxy (i.e., after they are 
accreted into a larger halo). 

4. THE PROPERTIES OF CENTRAL GALAXIES 

We now turn to a more detailed investigation of the 
properties of central galaxies in our group catalogue. 
We re-iterate that a central galaxy is either defined as 
the brightest group member (BCG) or the most massive 
group member (MCG). 

4.1. Central luminosity (stellar mass) - halo mass 
relation 

In Fig. [SI we show the luminosity - halo mass (left 
panel) and stellar mass - halo mass (right panel) rela- 
tions. The solid circles and shaded areas indicate the 
median and 68% confidence levels of these relations ob- 
tained using Ml as halo mass, while the open squares 
with error-bars are the results obtained using Ms as halo 
mass. Clearly, and not surprisingly, the luminosity (stel- 
lar mass) of the BCG (MCG) increases with halo mass. 
In both cases we find that the slope of the relation de- 
creases significantly with increasing halo mass, in good 
agreement with previous results (e.g.. Vale & Ostriker 
2004,2006; Cooray 2005; Yang et al. 2003; Y05c; van 
den Bosch et al. 2007a). The physical reason for this 
change in slope is probably a combination of AGN feed- 
back, and changes in the efficiencies of radiative cooling 
and dynamical friction (e.g. Lin et al. 2004; Dekel 2004; 
Cooray & Milosavljevic 2005). 

To quantify the Lc — Mh relation shown in the left 
panel, we fit the data using the following function, 



Lr — Lr 



[Mn/Mir+f' 
(1 + Mh/MiY 



(6) 



This model contains four free parameters: a normalized 
luminosity, L^, a characteristic halo mass. Mi, and two 
slopes, a and (3. The solid line shown in the left panel is 
the best fit to the average Lc — Mh relation. Note that 
using Ml or Ms as the halo mass does not lead to any 
significant changes in the results. The best fitting param- 
eters are [login, log Ml, a, [3] = [10.45, 12.54, 0.17, 0.51]. 
Thus, according to Eq. [6l Lc scales with Mh roughly as 
Lc oc M^-i"^ for halos with Mu > lO^^.s /j-IMq, and as 
Lc cx MO-68 for halos with Mh < lO^^-^ h'^UQ. Un- 
fortunately, since we do not have data for halos with 
Mh < 10^^'^ h~^MQ, any change of behavior at the low- 
mass end is not constrained. For the M*^c — Mh relation 
shown in the right panel, we use a similar function to fit 
the data: 

_ {Mh/M,r+f' 

^--*'"(1 + M,/Mi)/^- 

The solid line shown in the right panel is the best fit of 
this model to the data, and the best-fit parameters are 
[log Mo, logMi, a, 13] = [10.86, 12.08, 0.22, 1.61]. Thus, 
M*,c scales with Mh as M»_c cx M"-^^ for halos with 
Mh > lO^^-i /j-^Mq, and as M*,c oc M^-®^ ]^a,los 
with Mh < 10121 /i-IMq. 
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TABLE 1 

The best fit parameters of CLFs for all, red and blue galaxies 



Galaxy 


log[A4] 


log([A/h]) 


rs 




logic 


0"c 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



ALL 



14.40 
14.10 
13.80 
13.50 
13.20 
12.90 
12.60 
12.30 
12.00 



15.00) 
14.40) 
14.10) 
13.80) 
13.50) 
13.20) 
12.90) 
12.60) 
12.30) 



14.58 
14.23 
13.94 
13.64 
13.34 
13.05 
12.75 
12.45 
12.16 



35.51 ±3.88 
23.68 ± 1.07 
15.27 ±0.90 
9.60 ±0.82 
5.72 ±0.60 
3.27 ±0.60 
1.87 ±0.29 
1.12 ±0.23 
0.67 ±0.12 



-1.66 ± 0.11 
-1.44 ±0.04 
-1.33 ± 0.07 
-1.20 ±0.05 
-1.11 ± 0.04 
-1.08 ± 0.08 
-1.09 ± 0.06 
-1.07 ±0.09 
-1.10 ±0.06 



10.799 
10.714 
10.649 
10.584 
10.513 
10.442 
10.350 
10.224 
10.074 



±0.019 
±0.014 
±0.012 
± 0.006 
± 0.003 
±0.013 
±0.020 
±0.019 
±0.018 



0.141 
0.146 
0.157 
0.149 
0.144 
0.137 
0.128 
0.107 
0.107 



±0.021 
±0.011 
±0.005 
± 0.009 
± 0.009 
±0.012 
±0.012 
± 0.020 
± 0.026 



RED 



14.40 
14.10 
13.80 
13.50 
13.20 
12.90 
12.60 
12.30 
12.00 



15.00) 
14.40) 
14.10) 
13.80) 
13.50) 
13.20) 
12.90) 
12.60) 
12.30) 



14.58 
14.23 
13.94 
13.64 
13.34 
13.05 
12.75 
12.45 
12.16 



27.14 ±3.25 
18.72 ±2.53 
11.79 ±0.76 
6.98 ±0.80 
3.97 ±0.59 
2.12 ±0.44 
1.18 ±0.19 
0.67 ±0.15 
0.39 ±0.09 



-1.68 ± 0.13 
-1.43 ± 0.07 
-1.30 ±0.07 
-1.17 ±0.04 
-1.07 ±0.05 
-1.06 ± 0.09 
-1.03 ± 0.04 
-0.96 ± 0.09 
-0.84 ±0.16 



10.801 
10.709 
10.644 
10.581 
10.510 
10.438 
10.335 
10.207 
10.046 



±0.022 
±0.051 
±0.011 
± 0.005 
± 0.008 
±0.011 
±0.026 
±0.031 
± 0.046 



0.147 
0.144 
0.153 
0.146 
0.144 
0.140 
0.126 
0.100 
0.089 



±0.021 
±0.029 
±0.006 
± 0.008 
± 0.007 
±0.012 
±0.023 
±0.014 
± 0.005 



BLUE 



[14.40 
[14.10 
[13.80 
[13.50 
[13.20 
[12.90 
[12.60 
[12.30 
[12.00 



15.00) 
14.40) 
14.10) 
13.80) 
13.50) 
13.20) 
12.90) 
12.60) 
12.30) 



14.58 
14.23 
13.94 
13.64 
13.34 
13.05 
12.75 
12.45 
12.16 



9.21 ± 1.98 
4.94 ±0.57 
3.59 ±0.42 
2.51 ±0.27 
1.69 ±0.34 
1.08 ±0.24 
0.66 ±0.10 
0.43 ± 0.08 
0.28 ±0.06 



-1.55 ± 0.09 
-1.45 ± 0.05 
-1.37 ±0.08 
-1.26 ± 0.08 
-1.17±0.11 
-1.14 ±0.10 
-1.20 ±0.06 
-1.20 ±0.10 
-1.25 ± 0.09 



10.770 
10.740 
10.658 
10.603 
10.546 
10.473 
10.392 
10.277 
10.135 



± 0.043 
± 0.024 
± 0.024 
±0.017 
±0.019 
±0.016 
±0.022 
± 0.037 
±0.052 



0.130 
0.174 
0.191 
0.165 
0.166 
0.147 
0.131 
0.102 
0.098 



± 0.045 
±0.019 
±0.012 
±0.011 
± 0.009 
±0.011 
±0.015 
± 0.020 
±0.015 



Note. — Column (1): Galaxy sample. Column (2): halo mass range. Column (3): average of the 
logarithm of the halo mass. Column (4)-(7): average of the best fit free parameters to the four measurements 
of the CLFs, as shown in Fig. ^ The errors indicate the scatter among these four measurements or the 
scatter obtained from the 200 bootstrap samples, whichever is larger. 




12 13 14 
log[M,/(h-'Mj] 




12 13 14 
log[M,/(h-'M<,)] 



Fig. 6. — The left-hand panel shows the median luminosity of the BCG, Lc, as function of halo mass, Mf^. The solid dots and open 
squares show the results for halo mass and Ms, respectively. The related shaded areas (or error- bars) show the 68% confidence regions 
of Lc- The right-hand panel shows the median stellar mass of the MCG, M^^c, as function of halo mass, M^. Again solid dots and open 
squares show the results for halo mass Mj^ and Ms, and shaded area (or error-bars) show the 68% confidence regions of M^^c- The solid 
lines shown in the left and right panels are the best fit results using equations \E\ and [T] respectively. 
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12 13 14 15 12 13 14 15 

log[M,/(h->M3)] log[M,/(h-'M^)] 

Fig. 7. — The left-hand panel shows the luminosity gap between the BCG and the brightest satellite galaxy in the group, log Li — logL2i 
as function of group mass A/^. The right-hand panel shows the stellar mass gap between the MCG and the most massive satellite galaxy, 
logA/»_i — logAf«,2 as function of group mass Mf^. Results are shown for two cases, where groups with only 1 member are treated 
differently. In case 1, we assume that the brightest satellite galaxy has zero luminosity (open squares), while in case 2 we assume that the 
brightest satellite galaxy has an apparent magnitude equal to the magnitude limit of the survey (solid dots). The shaded areas indicate 
the corresponding 68% confidence intervals. The solid line in the left panel shows the best fit values of (T(log Lc) obtained from the CLFs. 
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Fig. 8. — The probability distributions of luminosity gap log Li — logL2 (upper panels) and stellar mass gap logM,,i — logM*_2 (lower 
panels) in halos of different mass ranges, as indicated. Similar to Fig. [Jl results are shown for both case 1 (unshaded histograms) and case 
2 (shaded histograms). In case 1, the isolated galaxies are put into the bin with logLi — logL2 = 2.0 or logM* i — log A/, 2 = 2.0. 
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4.2. The luminosity (stellar mass)-gap statistic 

An useful quantity to describe the difference be- 
tween the central and satellite galaxies is the so-called 
'luminosity-gap' between the BCG and the brightest 
satellite galaxy in a given halo. This 'luminosity-gap' 
holds important information regarding the formation 
and evolution of galaxies. For example, as discussed in 
D'Onghia et al. (2005) and Milosavljevic et al. (2006), 
the luminosity-gap may be used to quantify the dynam- 
ical age of a system of galaxies: halos with L2/L1 close 
to unity must be relatively young, as dynamical fric- 
tion causes multiple luminous galaxies in the same halo 
to merge on a relatively short time scale. Put differ- 
ently, the distribution of L2/L1 holds important infor- 
mation regarding the importance of galactic cannibalism 
for BCGs (e.g., Tremaine & Richstone 1977; Vale & Os- 
triker 2007) 

Using our group catalogue, it is straightforward to 
measure L2/L1, as well as the stellar mass equivalent, 
M^,2/M^*,i, as function of group mass, as long as the 
group has at least two members. In the case of 'iso- 
lated galaxies', i.e., groups with only one member (the 
central), however, there are two possibilities: either the 
galaxy is truly isolated, in that its halo contains no satel- 
lite galaxies, or the halo contains one or more satellite 
galaxies that are fainter than the flux limit of the survey. 
To bracket the extremes we proceed as follows. For iso- 
lated galaxies we either assume that the brightest satel- 
lite has zero luminosity (case 1) or that its luminosity is 
such that its apparent magnitude limit is equal to the 
magnitude limit of our sample (case 2). Fig. [7| shows 
the results thus obtained for both the luminosity gap 
(left-hand panel) and the stellar mass gap (right-hand 
panel). Symbols and shaded areas indicate the medians 
and 68% confidence intervals, respectively, and results 
are shown for both case 1 (open squares) and case 2 (solid 
dots). For groups with Mh > 10^^-^ /I'^M©, case 1 and 
2 yield identical results, simply because all these groups 
contain at least two member galaxies above the magni- 
tude limit of the survey. In the most massive halos, with 
Mh > IO'^^-^H-'^Mq, the median gap is logLi -logL2 ~ 
log Ml — logM2 ^0.2. The median luminosity or stel- 
lar mass gap can be reliably measured down to a halo 
masses of Alh ^ 10^^'^ h~^MQ, where the median values 
are logLi — logL2 ~ 0.3 and log Mi — logM2 ~ 0.4. For 
halos with Mh ^ 10^^'^ H^^Mq, however, cases 1 and 2 
yield very different results, indicating that the flux limit 
of the survey severely impedes our ability to accurately 
measure the gap statistics. Finally we emphasize that 
the luminosity-gap statistics obtained from Samples II 
and III separately are very similar, indicating that fiber 
collisions do not have a strong impact on the results pre- 
sented here. 

As pointed out by Tremaine & Richstone (1977; and 
references therein), the 'luminosity-gap' can be used to 
determine whether the BCGs in galaxy groups are sta- 
tistically drawn from the same distribution as the satel- 
lite galaxies or whether they are 'special'. If the av- 
erage magnitude difference between the BCGs and the 
brightest satellite galaxies is smaller than the standard 
deviation in the magnitudes of the BCGs, then they are 
consistent with being draw from the same distribution. 
To test this, the solid curve in the left-hand panel of 



Fig. [7] indicates the best-fit values of CTc = CT(logLc) ob- 
tained from the CLFs. A comparison with the median 
luminosity gaps suggests that the BCGs in groups, espe- 
cially with Mh ^ 10^^ °/i~^Mq, form a 'special' subclass, 
in that their luminosities can not be considered the ex- 
treme values of the distribution of satellite luminosities. 
As a cautionary remark, we emphasize that because of 
the method used to assign halo masses to the groups, 
the value of CTc may be underestimated, especially for 
low mass groups (see discussion in section [3]). However, 
even if CTc were underestimated by a factor two at the 
low mass end, this would not change our conclusion that 
BCGs are special in low mass halos. 

The 'luminosity-gap' can also be described using its 
distribution for groups of a given mass. In the upper 
panels of Fig. [8] we show the distribution of log L 1 — log L2 
for groups in the above mentioned case 1 (unshaded his- 
tograms) and case 2 (shaded histograms), respectively. 
Results are shown for three different mass bins, as indi- 
cated in the panels. In the lower panels of Fig.[8l results 
are shown for the corresponding stellar mass gap dis- 
tributions. One can see, the width of the distribution 
increases with decreasing halo mass. 

Systems with a relatively large luminosity gap, which 
most likely owes to the fact that the brightest galaxies in 
the halo have merged, are often termed "fossil groups" 
and have received a significant amount of attention in 
the recent literature (see Vikhlinin et al. 1999; Jones et 
al. 2003; D'Onghia et al. 2005; Milosavljevic et al. 2006; 
Sommer-Larsen 2006; van den Bosch et al. 2007a; Sales 
et al. 2007; von Benda-Beckmann et al. 2007). Follow- 
ing Jones et al. (2003) and Milosavljevic et al. (2006) 
we define systems in which the brightest satellite galaxy 
is at least 2 magnitudes fainter than the BCG (i.e., 
logLi — logi2 > 0.8, indicated as the dotted vertical 
lines in Fig. [8]), as "fossil" systems. From the SDSS 
DR4 group catalogue, we obtain that the fraction of fos- 
sil systems increases from '--^ 0.5 percent for groups with 
Alh ~ lO"^''-^ h~^MQ, to ~ 2.5 percent for groups with 
Mh 10"° /i-^Mq, 11 - 20 percent for groups with 
Mh ^ 10^^-^ /i"^Mq, and 18 - 60 percent for groups 
with Mh 10^^ ° /j-^Mq''. These resuhs are in good 
agreement with a similar analysis of galaxy groups in the 
2dFGRS by van den Bosch et al. (2007a). On the other 
hand, Jones et al. (2003) obtained an incidence rate of 
8 to 20 percent for systems with an X-ray luminosity 
from diffuse, hot gas of Lx,boi > 2.5 x lO^^ft-^^ergs^^. 
However, since the groups in our SDSS DR4 catalogue 
are not X-ray selected, a detailed comparison with our 
results is not possible. In a recent paper, D'Onghia et 
al. (2005) used detailed hydrodynamical simulations to 
predict the fraction of halos with Mh ^ lO^'*/i~^M0 that 
have logLi — logL2 > 0.8. From a total of twelve sim- 
ulated groups, they obtain a fossil fraction of 33 ± 16 
percent. This value is much higher than the fraction of 
fossil systems we find in the SDSS, which suggests a po- 
tential over-merging problem in their simulations. More 
recently von Benda-Beckmann et al. (2007) used a com- 
bination of N-body simulations and empirical models for 
the connection between galaxy luminosity and halo mass 
(taken from Cooray & Milosavljevic )2005), and found 

^ Whenever two values are quoted, these reflect the two extreme 
cases described above. 
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TABLE 2 

Halo occupation of satellite galaxies 



Satellite galaxies 


log A/s,o 


7 




(1) 


(2) 


(3) 


< 


-18.0 


12.48 ±0.04 


1.01 ±0.06 


< 


-18.5 


12.62 ±0.02 


1.05 ± 0.04 


< 


-19.0 


12.77 ±0.02 


1.06 ± 0.03 


< 


-19.5 


12.93 ±0.01 


1.07 ±0.02 


< 


-20.0 


13.15 ±0.01 


1.09 ± 0.02 


< 


-20.5 


13.44 ±0.01 


1.10 ±0.02 


< 


-21.0 


13.82 ±0.01 


1.13 ± 0.02 


< 


-21.5 


14.34 ±0.01 


1.33 ± 0.07 


Note. 


— Halo occupation model 


parameters 



for satellite galaxies in tlie SDSS DR4. Here the 
mean lialo occupations of satellite galaxies are de- 
scribed by (Ns) = (Mh/Ma,o)^. Column (1): The 
absolute magnitude limit of the satellite galax- 
ies. Columns (2-3): the best fit parameters Ms^o 
and 7 (averages with errors). The errors are esti- 
mated from the variances between Samples II and 
HI, with halo masses and Afg, and are much 
larger than the errors obtained from our 200 boot- 
strap samples. 

that the fossil group fraction is about 24% among all 
systems with masses 1 — 5 x lO^'^ ^ H^^Mq. This is in 
good agreement with our direct measurement from the 
SDSS data. 

5. HALO-OCCUPATION STATISTICS 

The upper panels of Fig. [9] show the mean halo oc- 
cupation numbers, (N) (the mean number of all group 
members) and (Ns) (the mean number of satellites), as 
functions of halo mass Mh- Results are shown only for 
Alfi — Ml, but adopting Mh = Ms gives very similar re- 
sults. The results shown in the left, middle and right pan- 
els correspond to galaxies in different luminosity ranges, 
as indicated in the panels. Similar to what is found in 
Y05c, the sharp break at the low mass end at {N) ~ 1 
^ indicates an almost deterministic relation between the 
luminosity of a central galaxy and the mass of its dark 
matter halo. The shoulder suggests that the brightest 
satellite galaxies is in general much fainter than the cen- 
tral galaxy (e.g. Zheng et al. 2005). The mean number 
of satellite galaxies, which is shown as the dashed line, 
reveals a very good power law feature. To quantify this, 
we model the halo occupation for satellite galaxies with 
{Ns) — {Mh/Msfi)'' , where Mg^ is a characteristic halo 
mass above which there is on average at least one satel- 
lite galaxies and 7 is the power law index. The best fit 
parameters for satellite galaxies with different absolute 
magnitude cuts are listed in Table [2] As an illustration, 
we show in the upper panels of Fig. [5] the corresponding 
best fit model predictions of the halo occupation of satel- 
lite galaxies as the dotted lines. They all agree with the 
data extremely well. 

In addition to the mean halo occupation number, we 
also investigate the second moment of the halo occupa- 
tion distribution (see Y05c). Here we only consider satel- 
lite galaxies, because by definition a central galaxy is al- 
ways assigned to a group in a deterministic way. This 
quantity is crucial in modelling the two-point correla- 

''' This is not seen in the upper left panel because our group 
catalogue does not include halos with masses below 10^^'® h Mq. 



tion function of galaxies on small scales (e.g., Benson et 
al. 2000; Berlind et al. 2003; Y05c; Tinker et al. 2007), 
and holds important information regarding the physical 
processes related to galaxy formation. In earlier inves- 
tigations, a number of simple models were adopted to 
describe the second moment of the halo occupation dis- 
tribution and its dependence on halo mass (e.g., Berlind 
& Weinberg 2002; Yang et al. 2003). In particular, us- 
ing the group samples constructed from the 2dFGRS, 
Y05c measured the second moment for all group mem- 
bers and found that the halo occupation distribution is 
close to Poissonian in massive halos and significantly sub- 
Poissonian in low mass halos. In the middle row of Fig. [9l 
we show {N'^) for satellite galaxies as a function of halo 
mass Mh- As one can see, {N^} is roughly proportional 
to M,j. To see how the distribution deviates from a Pois- 
son distribution, we show, in the lower panels of Fig [51 
{Ng) — (Ns)'^ as a function of (Ns). A Poisson distribu- 
tion has (N^) - {Ns)"^ = (Ns) (doted fines), while a deter- 
ministic distribution has {N^) — {Ngy ~ 0. The results 
thus indicate that the number distribution of satellite 
galaxies follows roughly a Poisson distribution. 

Fig.fTOlshows the distribution of the number of satellite 
galaxies in groups, Ng, for different halo mass bins. The 
thick solid curves indicate Poisson distributions with the 
same (Ns). As one can see, the observed A^s-distributions 
are well fitted by Poisson distributions. These properties, 
which have already been found in Y05c, suggest a direct 
link between satellite galaxies and dark matter sub-halos. 
In a recent study, Kravtsov et al. (2004), using large nu- 
merical simulations, have shown that the occupation dis- 
tribution of dark matter sub-halos follows Poisson statis- 
tics. This suggests that there may be a tight link between 
satellite galaxies and dark matter sub-halos, which has 
been a standard assumption in various HOD / CLF mod- 
els (e.g.. Vale & Ostriker 2004; 2006). Our resuffs provide 
observational support for such a link. 

6. SATELLITE FRACTIONS 

The satellite fraction as function of luminosity, Jsat{L), 
is an important quantity for a proper interpretation of 
measurements of galaxy-galaxy lensing (e.g., Guzik & 
Seljak 2002; Mandelbaum et al. 2006; Yang et al. 2006) 
and pairwise velocity dispersion of galaxies (e.g., Jing & 
Borner 2004; Slosar, Seljak & Tasitsiomi 2006; van den 
Bosch et al. 2007a; Li et al. 2007). In addition, since the 
halo bias depends on halo mass (Mo & White 1996), and 
since a satellite galaxy resides in a more massive halo 
than a central galaxy of the same luminosity (van den 
Bosch et al. 2007b), the large scale clustering of galax- 
ies of a given luminosity also depends strongly on the 
fraction of satellite galaxies. Here we estimate /sat(i) 
directly from our group catalogue. According to the test 
we carried out in the previous section, this fraction can 
be determined relatively accurately. 

In the left-hand panel of Fig. [11] we show fsat{L) as 
a function of galaxy luminosity. The results are plot- 
ted separately for all (solid lines) , red (dashed lines) and 
blue (dotted lines) galaxies. Since fiber collisions are ex- 
pected to significantly impact the number of close pairs, 
it can affect the satellite fractions fsntiL). To assess the 
uncertainty induced by the fiber collisions, we show re- 
sults for both Samples II and HI, as thick and thin lines, 
respectively. As mentioned earlier. Samples II and III 
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Fig. 9. — Upper row of panels: the average halo occupation numbers as function of halo mass for all galaxies {{N), solid lines) and 
for satellite galaxies only {{Ns), dashed lines). Here halo masses are taken to be Mi^. Panels on the left, in the middle and on the right 
correspond to an absolute mag nitudc limit of O ^Mr - Slog /i < -19.0, -20.0 and -21.0, respectively. The error-bars indicate the l-cr 
variances obtained from 200 bootstrap samples. Middle row of panels: similar to the upper row of panels, except that here we plot the 
second moment of the occupation numbers of satellite galaxies, {N^}. Lower row of panels: similar to the middle row of panels, except 
that now we plot (A''^) — (Ns)'^ as function of (Ng). The dotted, diagonal line indicates {N^) — (Ns)'^ = {Ns) and corresponds to Poisson 
statistics. Apparently, the distribution of Ns is very similar to a Poisson distribution. 



may respectively under- and over-estimate the number of 
group members because of their different treatments of 
fiber colhsions. Most likely the sateUite fractions of these 
two extreme cases bracket the true satehite fractions. As 
one can see, the satehite fraction decreases with increas- 
ing luminosity, from ^ 40% at "-^M^ - 5 log ft. = -17.0 
to ~ 5% at " iM^ - 51og/i = -22. The sateUite 
fraction of red galaxies at the faint end, ~ 70% at 
°-^Mr — 51og/i = —17.0, is significantly higher than that 
of blue galaxies, - 30% at °-'^Mr - 5 log /i = -17.0. The 
satellite fraction as a function of galaxy stellar mass is 
shown in the right panel of Fig. [TT] The overall behav- 
ior here is similar to that shown in the left panel. For 
reference, we list in Table [3] the fraction of satellites as a 
function of galaxy luminosity or stellar mass, separately 
for all, red and blue galaxies. The averages listed in this 
table are the averages between Sample II and III, while 



the errors quoted represent the deviations of the samples 
from the average. 

In recent years, the fraction of satellite galaxies 
has been studied extensively using HOD/CLF models 
(Cooray 2006; Tinker et al. 2007; van den Bosch et al. 
2007a), and analyses of galaxy-galaxy lensing measure- 
ments (e.g. Mandelbaum et al. 2006). For comparison, 
we overplot in Fig. [TT] the results obtained by Mandel- 
baum et al. for early-type galaxies (solid circles) and 
late-type galaxies (open triangles), with error-bars indi- 
cating the 95% confidence level. Although their samples 
are defined differently from ours (they separate galaxies 
into early- and late-types according to galaxy morpholo- 
gies, while we separate galaxies according to colour), our 
results match theirs quite well. 

7. SUMMARY 
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Fig. 10. — Number distributions of the satellite galaxies in groups of different halo mass bins, as indicated. Panels in the upper, middle 
and lower rows correspond to different absolute magnitude limits as indicated. The hatched histograms indicate the distributions obtained 
from the groups in the SDSS. Thick solid curves correspond to Poisson distributions with the same mean Ns, and are shown to illustrate 
the Poissonian nature of P{Ns\Mh)- 



Using a large galaxy group catalogue constructed from 
the SDSS Data Release 4 (DR4) by Y07, we have in- 
vestigated various halo occupation statistics of galaxies. 
In particular, we have split the galaxy population in red 
and blue galaxies, and in centrals and satellites, and de- 
termined the conditional luminosity functions of these 
varies subsamples. We have also presented luminosity 
gap statistics, satellite fractions, and halo occupation 
numbers for the galaxies in our group sample. The main 
results are summarized as follows: 

1. The conditional luminosity functions for central 
and satellite galaxies can be well modelled with a 
log-normal distribution and a modified Schechter 
form, respectively. The corresponding best fitting 
parameters are listed in Table 1. 

2. The average scatter of the log- normal luminos- 
ity distribution of central galaxies decreases from 



^ 0.15 dex at the massive end {log[Mh/ h ^Mq] > 
13.5) to ^ 0.1 dex at the low mass end 
(log[A'4/ H-^Mq] - 12.0). However, due to the 
method used to assign halo masses to the groups, 
at the low mass end this should be considered a 
lower limit on the true amount of scatter. 

3. The slope of the relation between the aver- 
age luminosity of central galaxies (in the "-^r- 
band) and halo mass, dlogLc/dlogM^,, decreases 
from ~ 0.68 for log[M?,/ /i-^M©] -C 12.5 to 
~ 0.17 forlog[Mfe/ /i-iM©] > 12.5. For the 
stellar masses of the central galaxies we ob- 
tain that dlog Af*^c/dlogAf;i, decreases from ^ 
1.83 for log[Mft/ h-^MQ] < 12.1 to ~ 0.22 
forlog[M,,/ /i-IMq] » 12.1. 

4. The halo (group) occupation numbers of satellite 
galaxies accurately follow Poisson statistics. Since 
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Fig. 11. — The fraction of satellite galaxies as function of luminosity (left-hand panel) and stellar mass (right-hand panel). Results are 
shown separately for all, red and blue galaxies as indicated. The thick lines are results for group Sample II and thin lines for group Sample 
III. For comparison, we show the satellite fractions obtained by Mandelbaum et al. (2006) from a galaxy-galaxy lensing analysis of the 
SDSS as the solid circles and open triangles with vertical errorbars (95% CL) for early- and late-type galaxies, respectively 



TABLE 3 
Satellite Fractions 



Galaxy 


All 


Red 


Blue 


(1) 


(2) 


(3) 


(4) 



-22.34 
-21.95 
-21.56 
-21.18 
-20.80 
-20.41 
-20.01 
-19.61 
-19.22 
-18.82 
-18.41 
-18.01 
-17.61 
-17.21 



11.28 
11.08 
10.89 
10.70 
10.50 
10.30 
10.10 
9.91 
9.70 
9.51 
9.31 
9.11 



0.021 
0.053 
0.100 
0.141 
0.186 
0.230 
0.264 
0.306 
0.336 
0.358 
0.372 
0.382 
0.396 
0.400 



± 0.001 
± 0.010 
± 0.012 
± 0.025 
± 0.040 
± 0.047 
± 0.052 
± 0.050 
± 0.056 
± 0.062 
± 0.054 
± 0.057 
± 0.064 
± 0.070 



0.023 
0.058 
0.116 
0.167 
0.230 
0.289 
0.340 
0.400 
0.453 
0.505 
0.548 
0.584 
0.624 
0.650 



± 0.001 
± 0.008 
± 0.013 
± 0.029 
± 0.045 
± 0.055 
± 0.062 
± 0.060 
± 0.067 
± 0.078 
± 0.069 
± 0.080 
± 0.083 
± 0.086 



0.017 
0.040 
0.069 
0.104 
0.135 
0.164 
0.185 
0.211 
0.233 
0.253 
0.269 
0.275 
0.280 
0.288 



± 0.003 



0.018 
0.011 
0.021 
0.033 
0.036 
0.039 
0.037 
0.040 
0.042 
0.033 
0.031 
0.035 
0.044 



0.061 
0.108 
0.155 
0.197 
0.233 
0.260 
0.282 
0.300 
0.304 
0.317 
0.328 
0.336 



± 0.009 
± 0.016 
± 0.029 
± 0.039 
± 0.047 
± 0.049 
± 0.052 
± 0.057 
± 0.057 
± 0.054 
± 0.054 
± 0.053 



0.062 
0.112 
0.164 
0.226 
0.287 
0.343 
0.401 
0.461 
0.507 
0.557 
0.587 
0.610 



± 0.008 
± 0.015 



0.030 
0.042 
0.051 
0.051 
0.052 
0.064 
0.055 
0.057 
± 0.066 
± 0.064 



0.049 
0.084 
0.123 
0.146 
0.166 
0.185 
0.199 
0.220 
0.237 
0.255 
0.266 
0.271 



± 0.026 



0.022 
0.022 
0.031 
0.041 
0.045 
0.046 
0.047 
0.050 
0.046 
0.041 
0.039 



Note. — Satellite fractions of galaxies in the SDSS DR4, 
as function of galaxy luminosity and stellar mass. Results are 
listed for all, red and blue galaxies, respectively. Column (1): 
The average luminosity (^'^M,, — 5 log h; upper part) and stel- 
lar mass {log[M^ /h~^ Mq]; lower part) of galaxies. Columns 
(2-4): the fractions of satellite galaxies for all, red and blue 
galaxies (averages with errors). The errors are estimated from 
the variances between Samples II and III, and are much larger 
than the errors obtained from our 200 bootstrap samples. 



the same applies to dark matter sub-halos, this sup- 
ports the standard picture that sateUite galaxies 
are associated with dark matter sub-halos. 



5. In massive halos with masses Mt > 10^"' h~^MQ 
roughly 85% (80%) of the central (satelhte) galax- 
ies are red. These red fractions decrease to 50% 
(40%) in halos with masses Mh 10^^ H-'^Mq. 

6. By comparing the scatter in the luminosities of 
BCGs to the luminosity difference between the 
EGG and its brightest satellite, we find that the 
BGGs form a 'special' subclass, in that their lumi- 
nosities can not be considered the extreme values of 
the distribution of satellite luminosities, expecially 
in halos with masses Mh < 10^^ ° H^'^Mq. 

7. The fractions of fossil groups, which are defined as 
groups a with luminosity gap logLi — logL2 > 0.8, 
decreases with increasing of halo mass from 18%- 
60% in halos with Mh lO^^ h~^MQ to ~ 2.5% in 
halos with Mh - lO^^^ /i^^Mq. 

8. The satellite fractions obtained from our group cat- 
alogue as functions of both luminosity and stellar 
mass (listed in Tablc[3|) are in good agreement with 
independent constraints from analyses of galaxy 
clustering and galaxy-galaxy lensing. 

These results can be used to constrain the various phys- 
ical processes related to galaxy formation and to inter- 
pret the various statistics used to describe large scale 
structures (e.g., galaxy correlation functions, pairwise 
velocity dispersions, etc.). Most of our findings are in 
good agreement with previous studies (e.g. Y05c, Zan- 
divarez et al. 2006; Robotham et al. 2006) and can be 
linked to the semi-analytical modelling of galaxy forma- 
tions (e.g., Kang et al. 2005; Zheng et al. 2005; Bower et 
al. 2006; Croton et al. 2006; De Lucia et al. 2007). The 
luminosity and stellar mass gap can be used to probe 
the specific formation properties of central galaxies (e.g., 
Vale & Ostriker 2007). The fraction of the red and blue 
populations for central and satellite galaxies can be used 
to probe the color evolution of satellite galaxies (ven den 
Bosch et al. 2007b). 



Galaxy Groups in SDSS DR4: II 



15 



ACKNOWLEDGMENTS 

We thank J. P. Ostrikcr and Zheng Zheng for helpful 
eomments, and Cheng Li for the galaxy color bi-normal 
fitting. XY is supported by the One Hundred Talents 
project, Shanghai Pujiang Program (No. 07pjl4102), 
973 Program (No. 2007CB815402), the CAS Knowledge 
Innovation Program (Grant No. KJCX2-YW-T05) and 
grants from NSFC (Nos. 10533030, 10673023). HJM 
would like to acknowledge the support of NSF AST- 
0607535, NASA AlSR-126270 and NSF IIS-0611948. 
Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating 
Institutions, the National Science Foundation, the U.S. 
Department of Energy, the National Aeronautics and 
Space Administration, the Japanese Monbukagakusho, 
the Max Planck Society, and the Higher Education 
funding Coimcil for England. The SDSS Web Site is 
y./ / www.sdss.org/ [ The SDSS is managed by the 



Astrophysical Research Consortium for the Participat- 
ing Institutions. The Participating Institutions are the 
American Museum of Natural History, Astrophysical In- 
stitute Potsdam, University of Basel, Cambridge Uni- 
versity, Case Western Reserve University, University of 
Chicago, Drexel University, Fermilab, the Institute for 
Advanced Study, the Japan Participation Group, Johns 
Hopkins University, the Joint Institute for Nuclear As- 
trophysics, the Kavli Institute for Particle Astrophysics 
and Cosmology, the Korean Scientist Group, the Chi- 
nese Academy of Sciences (LAMOST), Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astron- 
omy (MPIA), the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, Ohio State Univer- 
sity, University of Pittsburgh, University of Portsmouth, 
Princeton University, the United States Naval Observa- 
tory, and the University of Washington. 



REFERENCES 



Adelman-McCarthy J.K., et al., 2006, ApJS, 162, 38 

Baldry I.K., Glazebrook K., Brinkmann J., Ivezic Z., Lupton R.H., 

Nichol R.C., Szalay A.S., 2004, ApJ, 600, 681 
Baldry I.K., Balogh M.L., Bower R.G., Glazebrook K., Nichol R.C., 

Bamford S.P., Budavari T., 2006, MNRAS, 373, 469 
Bell E.F., Mcintosh D.H., Katz N., Weinberg M.D., 2003, ApJS, 

149, 289 

Benson A.J., Cole S., Frenk C.S., Baugh CM., Lacey C.G., 2000, 

MNRAS, 311, 793 
BerUnd A.A., Weinberg D.H., 2002, ApJ, 575, 587 
BerUnd A. A., Weinberg D.H., Benson A. J., Baugh CM., Cole S., 

Dave R., Frenk C.S., Jenkins A., Katz N., Lacey C.G., 2003, 

ApJ, 593, 1 
BerUnd A.A., et al. 2006, ApJS, 167, 1 

Berlind A. A., Kazin E., Blanton M.R.. Pueblas S., Scoccimarro R., 

Hogg D.H., 2007, preprint ( astro-ph/0610524 ) 
Blanton M. R. et al. , 2003a, Apj, 592, 8i9 
Blanton M. R. et al. , 2003b, AJ, 125, 2348 

Blanton M.R., Eisenstein D.J., Hogg D.W., Schlegel D.J., 

Brinkmann J., 2005a, ApJ, 629, 143 
Blanton M.R. et al. , 2005b, AJ, 129, 2562 

Bower R.G., Benson A.J.; Malbon R., Helly J.C, Frenk C.S., 
Baugh CM., Cole S., Lacey C.G., 2006, MNRAS, 370, 645 

Bullock J.S., Wechsler, R.H., Somerville R.S., 2002, MNRAS, 329, 
246 

Coil A.L., et al., 2006, ApJ, 638, 668 

CoUister A.A., Lahav O., 2005, MNRAS, 361, 415 

Cooray A., 2005, MNRAS, 364, 303 

Cooray A., 2006, MNRAS, 365, 842 

Cooray A., Milosavljevic M., 2005, ApJ, 627, L89 

Croton D.J., et al. 2006, MNRAS. 365. 11 

Dekel A., 2004, preprint ( astro-ph/0401503 ) 

De Lucia G., et al. 2007, MNnAa, 374, 809 

D'Onghia E., Sommer-Larsen J., Romeo A.D., Burkert A., 
Pcdersen K., Portinari L., Rasmussen J., 2005, ApJ, 630, L109 
Guzik J., Seljak U., 2002, MNRAS, 335, 311 

Hansen S.M., Sheldon E.S., Wechsler R.H., Kocster B.P., 2007, 

preprint ( arXiv:0710.3780 ) 
Jing Y.P., Mo H.J., Borner G., 1998, ApJ, 494, 1 
Jing Y.P., Borner G., Suto Y., 2002, ApJ, 564, 15 
Jing Y.P., Borner G., 2004, ApJ, 617, 782 

Jones L.R., Ponman T.J., Horton A., Babul A., Ebeling H., Burke 

D.J., 2003, MNRAS, 343, 627 
Kang X., Jing Y.P., Mo H.J., Borner G., 2002, MNRAS, 336, 892 
Kang X., Jing Y.P., Mo H.J., Borner G., 2005, ApJ, 631, 21 
Kauffmann G., White S.D.M., Guiderdoni B., 1993, MNRAS, 264, 

201 

Kravtsov A.V., Berlind A. A., Wechsler R.H., Klypin A. A., 
Gottlober S., AUgood B., Primack J.R., 2004, ApJ, 609, 35 

Li C, Kauffmann C, Jing Y.P., White S.D.M., Borner G., Cheng 
F.Z., 2006, MNRAS, 368, 21 

Li C, Jing Y.P., Kauffmann G., Borner G, Kang X., Wang L., 
2007, MNRAS, 376, 984 

Lin Y.T., Mohr J. J., 2004, ApJ, 617, 879 

Mcintosh D.H., Guo Y., Hertzberg J., Katz N., Mo H.J., van den 
Bosch E.G., Yang X., 2007, preprint (arXiv:0710.2157) 

Magliocchetti M., Porciani C, 2003, MNRAS, 346, 186 

Mandelbaum R., Seljak U., Kauffmann G., Hirata CM., 
Brinkmann J., 2006, MNRAS, 368, 715 

Marinoni C, Hudson M.J., 2002, ApJ, 569, 101 



Milosavljevic M., Miller C.J., Furlanetto S.R., Cooray A., 2006, 

ApJ, 637, L9 
Mo H.J., White S.D.M., 1996, MNRAS, 282, 347 
Mo H.J., Yang X.H., van den Bosch F.C, Jing Y.P., 2004, MNRAS, 

349, 205 

Peacock J. A., Smith R.E., 2000, MNRAS, 318, 1144 
Petrosian V., 1976, ApJ, 209, LI 

Robotham A., Wallace C, Phillipps S., De Propris R., 2006, ApJ, 
652, 1077 

Sales L.V., Navarro J.F., Lambas D.G., White S.D.M., Croton D.J., 

2007, preprint (arXiv0706.2009) 
Schlegel D.J., Finkbeiner D.P., Davis M., 1998, ApJ, 500, 525 
Scoccimarro R., Sheth R.K., Hui L., Jain B., 2001, ApJ, 546, 20 
Scranton R., 2002a, MNRAS, 332, 697 
Seljak U.,2000, MNRAS, 318, 203 

Slosar A., Seljak U., Tasitsiomi A., 2006, MNRAS, 366, 1455 
Sommer-Larsen J., 2006, MNRAS, 369, 958 
Strauss M.A., et al., 2002, AJ, 124, 1810 

Tinker J.L., Weinberg D.H., Zheng Z., Zehavi I., 2005, ApJ, 631, 
41 

Tinker J.L., Norberg P., Weinberg D.H., Warren M.S., 2007, ApJ, 
659, 877 

Tremaine S.D., Richstone D.O., 1977, ApJ, 212, 311 

Vale A., Ostriker J. P., 2004, MNRAS, 353, 189 

Vale A., Ostriker J. P., 2006, MNRAS, 371. 1173 

Vale A., Ostriker J. P., 2007, preprint (astro-ph/0701096 ) 

van den Bosch F.C, Yang X., Mo H.J., 2003, MNnAa, 340, 771 

van den Bosch F.C, Mo H.J., Yang X., 2003, MNRAS, 345, 923 

van den Bosch F.C, Weinmann S.M., Yang X., Mo H.J., Li C, 

Jing Y.P., 2005, MNRAS, 361, 1203 
van den Bosch F.C, Yang X., Mo H.J., Weinmann S.M., Maccio 

A., More S., Cacciato M., Skibba R., Kang X., 2007a, MNRAS, 

376, 841 

van den Bosch F.C, Aquino D., Yang X., Mo H.J., Pasquali 

A., Mcintosh D.H., Weinmann S.M., Kang X., 2007b, preprint 

(arXiv:0710.3164) 
Vikhilin A., McNamara B.R., Hornstrup A., Quintana H., Forman 

W., Jones, C, Way C, 1999, ApJ, 520, LI 
von Benda-Beckmann A.M., D'Onghia E., Gottloeber S., Hoeft 

M., Khalatyan A., Klypin A., Mueller V., 2007, preprint 

(arXiv:0710.1297) 
Weinmann S.M., van den Bosch F.C, Yang X., Mo H.J., 2006a, 

MNRAS, 366, 2 

Weinmann S.M., van den Bosch F.C, Yang X., Mo H.J., Croton 

D.J., Moore, B., 2006b, MNRAS, 372, 1161 
Wang Y., Yang X.H., Mo H.J., van den Bosch F.C, Chu Y., 2004, 

MNRAS, 353, 287 
White M., Zheng Z., Brown M.J.I., Dey A., Jannuzi B.T., 2007, 

ApJ, 655, 69 

Yan R., Madgwick D.S., White M., 2003, ApJ, 598, 848 

Yan R., White M., Coil A.L., 2004, ApJ, 607, 739 

Yang X., Mo H.J., van den Bosch E.G., 2003, MNRAS, 339, 1057 

Yang X., Mo H.J., Jing Y.P., van den Bosch E.G., Chu Y.Q., 2004, 

MNRAS, 350, 1153 
Yang X., Mo H.J., van den Bosch E.G., Jing Y.P., 2005a, MNRAS, 

356, 1293 

Yang X., Mo H.J., van den Bosch E.G., Jing Y.P., 2005b, MNRAS, 

357, 608 

Yang X., Mo H.J., Jing Y.P., van den Bosch E.G. 2005c, MNRAS, 

358, 217 {Y05c) 



16 



Yang et al. 



Yang X., Mo H.J., van den Bosch F.C., Wcinmann S.M., Li C, 

Jing Y.P., 2005d, MNRAS, 362, 711 
Yang X., Mo H.J., van den Bosch F.C., 2006, ApJ, 638L, 55 
Yang X., Mo H.J., van den Bosch F.C., Pasquali A., Li C, Bardcn 

M., 2007, ApJ, 671, 153 (paper 1, Y07) 
Zandivarcz A., Martnez H.J., Merchan M.E., 2006, ApJ, 650, 137 
Zehavi L, ct al., 2002, ApJ, 571, 172 
Zehavi L, et al., 2004 ApJ, 608, 16 



Zehavi L, et al., 2005 ApJ. 630. 1 

Zheng Z., Tinker J.L., Weinberg D.H., Berlind A.A., 2002, ApJ, 
575, 617 

Zheng Z., et al., 2005, ApJ, 633, 791 

Zheng Z., Coil A.L., Zehavi I., 2007, ApJ, 667, 760 

Zheng Z., Weinberg D.H., 2007, ApJ, 659, 1 



